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RESEARCH  SUMMARY 

This  paper  reformulates  the  empirical  fire  spread 
equations  of  Rothermel,  based  on  a  tenfold  increase  in 
the  data  base  of  laboratory  test  fires.  The  experiment 
was  performed  at  zero  windspeed,  with  zero  slope  and 
uniform  fuel  particle  size.  Major  revisions  include  the 
propagating  flux  ratio,  the  reaction  velocity,  and  the 
moisture  damping  coefficient.  These  parameters  re- 
main functions  of  the  observable  fuel  bed  variables. 
Also,  we  redefine  or  reinterpret  many  of  the  fuels'  de- 
scriptive parameters  such  as  the  net  loading,  the  heat 
of  combustion,  fuel  bed  bulk  density,  flame  residence 
time,  and  heat  of  ignition.  As  before,  spread  rate  is 
based  on  an  energy  balance  ratio  of  power  source 
and  heat  sink  terms. 

We  use  Frandsen's  (1973)  formulation  for  effective 
bulk  density  to  account  for  that  fraction  of  the  fuel  bed 
that  must  be  raised  to  ignition  temperature.  The  en- 
thalpic  load  of  ignition  also  includes  the  energy  re- 
quired for  complete  pyrolyzation  of  the  fuel  as  reported 
earlier. 

We  also  use  Susott's  (1984)  method  of  determining 
the  "flaming"  heat  of  combustion,  thus  separating  the 


gaseous  and  carbonaceous  mass  fractions  of  pyroly- 
sis.  This  provides  that  the  heat  of  combustion  will 
account  for  many  of  the  physiochemical  effects  of  min- 
erals, and  so  on,  on  the  production  rate  of  the  gaseous 
fuel.  Thus  the  reaction  intensity  is  of  the  flames  alone 
and  specifically  excludes  energy  derived  from  burning 
char  whether  in  the  flaming  zone  or  not.  It  is  still 
measured  as  an  energy  density  within  the  combustion 
zone. 

In  this  new  formulation,  moisture  damping  is  com- 
pletely divorced  from  fire  extinction.  Fire  extinction 
is  properly  handled  by  a  probability  function  acting 
on  the  reaction  velocity.  Moisture  damping  is  a  simple, 
negative  exponential  function  of  the  fuel  moisture, 
wherein  each  fuel  sample  has  a  characteristic  mois- 
ture content.  In  the  experiment,  the  characteristic  fuel 
moisture  was  dependent  only  on  the  size  of  the  fuel 
particles;  however,  it  is  further  shown  that  it  will  vary 
greatly  with  the  physiology  of  different  fuels  (bark  cov- 
ered, live  or  fresh-fallen  foliage,  decayed  wood,  and 
so  on). 

The  formula  for  reaction  velocity  begins  with 
Anderson's  8d  (1985)  rule-of -thumb,  which  provides 
that  the  burning  residence  time  in  minutes  is  equal  to 
eight  times  the  fuel  particle  thickness  in  inches.  The 
burning  rate  is  limited  by  the  free  path  for  radiation  or 
heat  flux  to  propagate  through  the  fuel  bed,  and  further 
restricted  by  limiting  the  effect  of  excessive  fuel  load- 
ing or  bed  depth.  A  probability  function  shuts  the  fire 
down  at  very  light  fuel  loads  or  high  fuel  moisture  con- 
tent as  previously  reported.  The  concept  of  optimum 
packing  ratio  is  abandoned. 

Following  Rothermel,  the  propagating  flux  is  some 
fraction  of  the  reaction  intensity,  empirically  formulated 
as  the  propagating  flux  ratio.  It  is  an  elementary  func- 
tion of  the  free  optical  path  through  the  fuel  bed. 

Procedure  is  outlined  whereby  the  manager  in  the 
field  may  estimate  the  probabilities  that  his  real  fire 
behavior  will  depart  from  the  prediction. 
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INTRODUCTION 

Fire  remains  a  complex,  little-understood  phe- 
nomenon; however,  wildfire  behavior  prediction 
is  an  increasingly  important  requirement  of  land 
managers.  Fire  intelligence  is  needed  in  a  broad 
range  of  applications,  from  real-time  fire  behavior 
predictions  for  immediate  tactical  use,  to  long-range 
planning,  budgeting,  and  resource  allocation.  For- 
est and  range  fires  are  extremely  complex  interac- 
tions of  physical,  chemical,  and  thermodynamic 
processes  for  which  satisfactory  scientific  theories 
do  not  exist  or  are  far  from  complete.  Fire  research- 
ers and  model  builders  have  therefore  made  heroic 
efforts  to  develop  semiempirical  mathematical  mod- 
els that  anticipate  management  needs.  Even  then 
the  mathematical  models  are  commonly  stretched 
beyond  their  capabilities  by  enterprising  managers 
who  will  improvise  tools  that  might  solve  an  imme- 
diate problem. 

Because  a  precise  theory  is  not  available  and 
because  the  range  of  application  is  increasing  in 
breadth  and  the  information  becoming  more  impor- 
tant and  valuable,  empirical  mathematical  fire  mod- 
els must  be  as  dependable,  reliable,  accurate,  and 
robust  as  the  current  state  of  technology  can  make 
them.  This  is  the  purpose  of  the  research  reported 
here — to  enhance  the  versatility  and  accuracy  of 
existing  fire  behavior  models.  Basically,  this  is  a 
reworking  of  the  Rothermel  (1972)  fire  model  equa- 
tions. The  primary  objective  of  the  experiment  was 
to  reexamine  the  fuel  moisture  damping  effects  and 
the  marginal  limits  of  combustion  near  extinction. 
The  reformulation  of  the  fire  spread  equations  here 
should  be  viewed  as  a  secondary  result  of  the  experi- 
ment where  the  new  equations  become  hypotheses 
that  require  further  testing. 

The  equations  introduced  here  are  based  on  a 
tenfold  increase  in  the  data  base  of  laboratory  test 
fires.  This  experiment  was  expressly  designed  to 
cover  the  range  of  all  the  physical  (geometric)  vari- 
ables, fuel  moisture  content,  and  fuel  loadings  that 
might  be  encountered  in  surface  fuel  beds  in  the 
field.  The  revised  equations  are  expressed  as 


functions  of  the  fuel  bed  packing  ratio,  (3;  the  fuel 
particle  surface  area-to- volume  ratio,  a;  the  fuel 
bed  depth,  8;  and  the  fuel  moisture  content,  Mr 
The  experiment  was  performed  at  zero  windspeed, 
with  zero  slope,  and  uniform  fuel  particle  size.  Fuel 
particles  ranged  from  0.06  cm  to  1.27  cm  in  thick- 
ness; however,  each  fuel  bed  was  constructed  of  a 
single,  unmixed  size  class.  The  tabulated  experi- 
mental data  are  accessible  from  computer  files  at 
the  Intermountain  Fire  Sciences  Laboratory  (IFSL), 
Missoula. 

We  used  heuristic  reasoning  and  statistical 
"because  it  works"  types  of  argument  to  reach 
"conclusions"  and  to  generate  empirical  equations 
to  describe  events.  For  example,  our  procedure  for 
forming  the  expectation  value  for  spread  rate  as  the 
product  of  other  primary  "expectation  values"  is  not 
vigorously  legitimate  because  those  expectation 
values  were  not  derived  independently.  Thus,  this 
exposition  cannot  be  labeled  as  good  physical  the- 
ory; knowledge  of  wildfire  phenomena  has  not  yet 
reached  that  stage. 

New  empirical  formulas  were  constructed 
through  judicious  use  of  modified  multiple  linear 
regression — analysis  of  variance  techniques  with 
constraints  and  transformation  of  variables  imposed 
to  conform  to  preconceived  notions  of  physical  prin- 
ciples. We  examined  most  possible  combinations  of 
the  independent  fuel  bed  variables,  but  constrained 
the  functions  of  these  variables  to  "reasonable"  for- 
mulations that  might  be  expected  from  any  future 
physical  theory  or  that  could  be  easily  revised  by 
further  analysis.  That  is,  the  rigorous  statistical 
curve-to-data  fitting  is  here  tempered  by  arbitrarily 
limiting  the  independent  variables  to  those  that 
"ought  to  have"  a  proper  physical  effect  and  by  in- 
corporating much  of  the  extensive  practical  knowl- 
edge of  several  colleagues  whose  experience  with 
field  applications  has  given  guidance  to  how  fire 
models  should  behave.  The  chosen  combinations 
of  independent  variables  are  about  the  most  effi- 
cient descriptors  of  the  dependent  fire  parameters 
that  make  good  physical  sense. 
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Similarly,  most  of  the  numerical  coefficients  are 
truncated  to  show  only  meaningful  digits.  In  some 
cases,  the  experimental  statistics  indicated  that 
more  significant  digits  were  appropriate,  but  the 
practical  physics  overruled,  for  example,  an  expo- 
nent was  evaluated  to  be  0.486  ±  0.011,  but  in  the 
formulation,  0.5  (square  root)  makes  more  sense. 
The  data  would  confirm  that  the  two  numbers  were 
not  really  different  at  some  reasonable  significance 
level;  so,  the  square  root  is  used  for  practical  rea- 
sons! Note  that  this  problem  would  be  superfluous 
if  we  had  theoretically  hypothesized  a  square  root 
function  before  the  experiment  because  the  data 
would  have  confirmed  it.  We  are  in  the  philosophi- 
cal position  of  using  this  publication  to  form  a  new 
hypothesized  fire  spread  model  that  must  be  con- 
firmed or  rejected  by  further  testing  and  validation. 
Prior  wildland  fire  models  were  "confirmed"  by  sat- 
isfactory performance  in  field  application,  and  that 
procedure  is  anticipated  here. 

We  have  used  metric  measure  exclusively,  and 
dimensionless  ratios  are  expressed  as  fractions. 

DISCUSSION  OF  EXPERIMENTAL 
TEST  METHODS  AND  RESULTS 

Rothermel's  model  is  the  most  widely  applicable 
and  robust  of  the  current  wildland  fire  spread 
models.  The  reader  should  be  familiar  with  the 
Rothermel  formulation.  We  have  found  no  reason 
or  justification  in  the  extended  data  set  to  change 
Rothermel's  basic  formulation,  in  which  spread  rate 
is  based  on  an  energy  balance  ratio  of  power  source 
and  heat  sink  terms: 

Power  density  of  propagating  flux 

R  =  

Preheating  energy  of  fuel 

IP  (1) 

=  ~p^7 

where  Ipis  the  fire  propagating  intensity  provided 
by  the  heat  source,  QT  is  the  heat  sink  energy  re- 
quired to  heat,  dry,  and  pyrolyze  the  fuel,  and  p6 
is  the  fuel  bed  bulk  density. 

We  reformulate  the  empirical  equations  by  fitting 
the  curves  to  the  expanded  data  set  including:  the 
propagating  flux  ratio    the  reaction  velocity,  r,  and 
the  moisture-damping  coefficient,  T|M.  These  para- 
meters remain  functions  of  the  observable  fuel  bed 
variables.  Also,  we  redefine  or  reinterpret  many  of 
the  fuels'  descriptive  parameters  such  as  wn, 
x  (residence  time),  and  QT,  which  plays  an  equiva- 
lent role  to  the  formerly  used  Qi  . 


Variable  List 

We  list  here  the  complete  set  of  variables  that  are 
used  in  the  revised  fire  spread  equations. 

The  fuel  bed  descriptors  include: 

1.  Fuel  ovendry  loading,  wo,  kg/m2 

2.  Fuel  particle  size:  ratio  of  surface  area-to- 
volume,  a,  cm-1 

3.  Fuel  bed  depth,  8,  cm 

4.  Fuel  moisture  content,  Mf,  fraction  ovendry 
weight 

5.  Wind  speed,  U,  m/min 

6.  Slope,  tan(<))),  were  invariant,  zero,  for  this  ex- 
periment (for  completeness,  one  should  in- 
clude the  ambient  air  temperature,  but  its 
effect  is  very  small  and  thus  assumed  to  be 
constant;  25  or  30  °C  is  adequate;  similarly, 
for  relative  humidity.  These  parameters  have 
a  significant  impact  in  models  for  the  diurnal 
variation  of  fuel  moisture.  Here,  however, 
fuel  moisture  is  input  directly  and  is  not  par- 
ameterized). 

The  measured  thermochemical  properties  of  the 
fuels  include: 

7.  Heat  of  combustion  of  pyrolyzate  gases,  hv, 
kJ/kg  (hu  is  measured  as  the  combustion  en- 
ergy of  the  pyrolyzate  gases  per  initial  unit 
mass  of  the  solid  dry  fuel.  The  total  heat  of 
combustion,  hQ2,  kJ/kg,  was  measured,  com- 
pared, and  rejected  in  the  analysis  of  variance 
procedure).  The  mineral  content  was  meas- 
ured, but  found  unnecessary  when  the  burn- 
ing rates  and  intensities  are  based  on  the 
heat  of  combustion,  hv,  of  the  vaporized  frac- 
tion of  fuel. 

8.  Fuel  particle  ovendry  density,  pp,  kg/m3 

9.  Heat  for  pyrolysis,     ,  kJ/kg 

10.   A  characteristic  fuel  moisture  content,  M  , 
fraction  ovendry  weight,  is  used  to  character- 
ize moisture  damping,  replacing  the  moisture 
of  extinction,  M  ,  used  in  Rothermel's  1972 
model. 

In  addition,  we  have  a  few  new  numeric  parame- 
ters in  the  equations  below  that  may  vary  with  the 
thermochemical  character  of  a  different  fuel  type. 
Specific  comment  is  made  where  appropriate. 

The  foregoing  independent  variables  are  also  the 
input  variables  for  the  revised  spread  equations  and 
are  intended  to  be  comparable  to  Rothermel's  (1972, 
page  27)  list  of  input  parameters  that  definitively 
describe  the  fuel  bed  and  burning  conditions.  For 
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convenience,  the  following  parameters  are  explicitly 
defined: 

11.  Fuel  bed  bulk  density,  p6  =  wolh,  is  the  oven- 
dry  mass  per  unit  volume  of  fuel  bed;  concept 
and  definition  is  unchanged 

12.  Fuel  bed  packing  ratio,  p*  =  p6  /  pp, 
unchanged 

The  measured  dependent  variables  in  the  ex- 
periment are,  of  course,  the  primary  parameters 
for  which  the  predictive  model  is  needed.  These 
included: 

13.  Rate  of  spread,  R,  m/min 

14.  The  reaction  intensity,  IR,  kw/m2,  is  experi- 
mentally calculated  by  the  product  of  the 
heat  of  combustion  and  weight  loss  rate. 
The  mass  loss  of  the  combustion  zone  was 
measured  continuously  during  each  fire. 
Then  the  time  average  weight  loss  rate  per 
unit  area  of  combustion  zone  was  computed 

15.  Depth,  dp  cm,  of  flaming  zone  in  the  longitu- 
dinal direction  of  R  (to  determine  area  of 
flaming  zone  and  the  reaction  time  which 

is  critical  to  reaction  velocity) 

16.  Flame  length,  lp  cm,  for  correlation  checks 
on  burning  rates  and  intensity  parameters. 
Also,  fuel  modelers  need  to  predict  scorch 
heights,  model  crown  fires,  etc. 

Heat  Load 

Frandsen  (1973)  rationalized  that  the  entire  mass 
of  large  fuel  particles  did  not  rise  to  ignition  tem- 
perature at  the  onset  of  flaming  combustion;  that 
it  is  sufficient  to  pyrolyze  the  surface  only  deep 
enough  to  produce  pyrolyzate  gas  fuel  to  sustain 
flame.  The  effect  of  heating  only  part  of  the  fuel 
mass  was  interpreted  in  terms  of  effective  bulk 
density,  which  was  defined  as  that  fraction  of  the 
total  fuel  bed  bulk  density  that  must  be  raised  to 
ignition  temperature.  Thus,  Frandsen's  effective 
heating  number,  e,  is  defined  as: 

pb,  effective 

p6,  total 

and  experimentally  found  by  him  to  be 

<e>=  exp(-l/(0.22 ■  a))  (2) 
for  a  measured  in  cm-1. 

The  net  fuel  mass  that  must  be  heated  to  ignition 
is  equivalent  to  a  thin,  2.2-mm,  surface  layer  on 
the  thermally  thick  fuel  particles,  according  to 
Frandsen. 


The  onset  of  flaming  combustion  occurs  when  the 
woody  forest  fuels  reach  about  320  °C.  Even  though 
Frandsen  has  shown  that  the  entire  fuel  particle  is 
not  raised  to  ignition  temperature  uniformly,  but 
only  the  surface  of  the  fuel  particle  must  be  ignited, 
we  have  found  that  the  ignition  enthalpic  load  of  the 
fuel  mass  must  include  the  energy  required  to  com- 
plete pyrolyzation  of  the  fuel.  Susott  (1982b,  1984) 
has  shown  that  pyrolyzation,  that  is,  flammable  gas 
production,  is  essentially  complete  very  near  400  °C. 
The  energy  of  fuel  ignition,  QT ,  is  newly  defined  as: 

QT  =  Qf+MfQM  (3) 

where 

T400 

Qf=)  (dQ/dDdT 
amb 

and  the  heat  to  vaporize  the  fuel  moisture  is 

QM  =  4. 18(100  -  Tmb  +  540),  kJ/kg 

and  where  Susott's  procedure  provides  values  of 
dQ  I  dT  to  400  °C  and  above,  where  Dunlap's  spe- 
cific heat  equation  (as  used  by  Frandsen  and  Roth- 
ermel)  fails  dramatically.  Note  that  we  have  just 
circumvented  a  conflict  in  semantics;  QT  is  no  longer 
the  heat  load  of  "preignition"  and  the  term  "heat  of 
pyrolysis"  has  explicit  meaning  elsewhere;  so  we 
will  use  "heat  of  ignition"  for  the  total  heat  load, 
QT.  We  hope  this  terminology  will  be  less  confusing 
to  later  users. 

Reaction  Intensity 

The  concept  of  propagating  flux  that  drives 
the  fire  is  significantly  changed  from  Rothermel's 
original  description.  We  propose  that  the  energy 
released  from  only  the  flaming  area  and  due  to 
burning  of  the  pyrolyzate  gases  is  the  major  source 
of  power  density  that  drives  the  fire.  The  propagat- 
ing intensity  is  still  calculated  as  a  fraction,  of 
reaction  intensity:  Ip  =  K'  Ir>  but  here  IR  is  really  the 
reaction  intensity  of  the  flames  and  specifically  ex- 
cludes energy  derived  from  burning  char  whether 
in  the  flaming  zone  or  not.  It  is  still  measured  as 
an  energy  release  rate  per  unit  area  within  the 
combustion  zone. 

Susott  (1982a)  and  Wilson  (1985,  1987)  argue  that 
the  flaming  pyrolyzate  gases  are  the  primary  energy 
source  for  propagation  of  the  fire  front,  and  Susott 
has  developed  techniques  to  measure  the  gaseous 
heat  of  combustion,  hv.  Because  Susott's  method 
determines  the  "flaming"  heat  of  combustion  per 
unit  of  total  (initial)  fuel  mass,  thus  separating  the 
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gaseous  and  carbonaceous  fractions,  we  have  the 
additional  advantage  that  hu  accounts  for  many 
of  the  physiochemical  effects  (minerals,  etc.)  on 
the  mass  fraction  and  on  the  production  rate  of 
the  gaseous  fuel.  For  example,  the  effects  of  miner- 
als that  retard  pyrolysis  to  form  less  volatiles  and 
more  char  are  reflected  in  the  lower  measured  hv. 
Although  they  vary  widely  in  woody  fuels,  the  py- 
rolysis gases  account  for  about  60  percent  of  the 
total  combustion  energy  and  about  80  or  90  percent 
of  the  total  fuel  mass  loss.  While  we  have  settled 
on  the  gaseous  heat,  hv,  for  use  in  modeling  the  re- 
action intensity  for  the  fire  spread  equations,  this 
does  not  preclude  using  the  total  heat  of  combust- 
ion, h02,  for  many  other  applications.  For  example, 
it  may  be  appropriate  in  "fire  effects"  work  to  use 
h02  to  calculate  the  total  heat  released  to  the  site 
which  includes  the  burnout  of  the  residual  char 
fraction.  Thus  the  reaction  intensity  is  defined: 

IR  =  -hv  (dwjdt).  (4) 

The  other  terms  in  the  reaction  intensity  that 
must  be  experimentally  measured  are  the  mass 
loss  rate  per  unit  area  of  the  combustion  zone  and 
the  reaction  or  residence  time.  Here,  again,  our 
method  differs  from  that  of  Rothermel  (1972).  In 
this  experiment,  we  measured  the  mass  loss  rate 
of  the  entire  combustion  zone  (ho,  is  measured  rela- 
tive to  total  mass)  and  the  width  and  depth  of  the 
burning — flaming  area.  Thus  a  consistent  new  cal- 
culation of  flaming  reaction  intensity  follows  as  the 
average  heat  release  rate  per  unit  area  of  the  flam- 
ing zone.  The  (flaming)  reaction  time  is  calculated 
directly  by  dividing  the  flame  zone  depth  by  the 
spread  rate. 

Note  that  Rothermel's  formulation  is  consistent 
in  the  same  sense:  he  directly  measured  the  mass 
loss  rate  per  unit  area  of  the  combustion  zone,  that 
when  multiplied  by  the  total  heat  of  combustion, 
gives  a  reaction  intensity  equal  to  the  total  heat 
release  rate  per  unit  area  of  combustion  zone.  The 
reaction  time  (and  hence  the  combustion  zone  area) 
then  was  determined  by  the  duration  of  constant 
mass  loss  in  the  combustion  zone.  The  difference 
between  the  two  methods  is  that  Rothermel's  com- 
bustion zone  was  defined  by  the  area  of  constant 
mass  loss,  while  here  it  is  the  area  with  attached 
flame. 

In  order  for  the  experimentalist  to  express  the 
reaction  intensity  as  a  function  of  the  fuel  bed  de- 
scriptors, it  is  convenient  to  make  three  simplifying 
assumptions: 

1.  The  reaction  intensity  is  linearly  proportional 
to  the  heat  of  combustion  of  the  fuel  and  to  the  loss 


rate  of  fuel  mass  per  unit  area,  as  in  equation  4. 
(Equation  4  actually  defines  the  reaction  intensity.) 
This  makes  the  reaction  intensity  dimensionally 
equivalent  to  a  power  flux  density  rising  from 
the  combustion  zone.  We  use  an  independently 
determined  hv,  which  incorporates  any  effect  of 
incomplete  combustion  of  the  pyrolyzates  and  other 
variations  in  hv  caused  by  the  nature  of  the  fire 
or  other  physical  parameters. 

2.  The  loss  rate  of  fuel  loading  is  linearly  propor- 
tional to  the  fuel  load.  The  "constant"  of  proportion- 
ality (that  Rothermel  dubbed  the  "reaction  velocity, 
r")  is  experimentally  fitted  to  a  function  of  the  ob- 
servable fuel  bed  characteristics.  Note,  as  shown 
below,  that  this  reaction  velocity  has  a  residual 
dependence  on  fuel  loading,  which  means  that  this 
assumption  that  the  fuel  loss  is  linearly  propor- 
tional to  fuel  loading  is  not  precisely  correct. 

IR  =  -hu(dwo/dt)  =  -hvwor  (5) 

Recall  that  Rothermel  used  the  net  loading,  wn,  to 
provide  for  less  than  100  percent  burning  efficiency, 
and  which  fraction  of  the  total  load,  w  ,  is  to  be  ac- 
counted  for  in  the  reaction  velocity,  T,  along  with 
the  reaction  time,  x.  In  the  present  experiment,  the 
residual  load  after  complete  burnout  never  exceeded 
0.005  fraction  of  the  initial  load,  so  by  Rothermel's 
method,  the  net  and  total  loads  were  equivalent  in 
the  experiment.  Also,  in  field  application,  fuel  mod- 
elers are  well  versed  in  selecting  "that  portion  of  the 
fuel  bed  that  will  contribute  to  propagating  the  fire." 
Our  concept  here  is  that  all  inefficiencies  in  flame 
production  will  be  accounted  for  by  Susott's  tech- 
nique for  hv  and  that  the  reaction  velocity  should 
concern  only  the  time  dependence  of  the  burning 
rate  function.  We  wish  to  be  explicit  that  the  carbo- 
naceous char  in  many  fires  was  almost  completely 
consumed  in  the  flaming  zone,  thus  contributing 
to  the  total  mass  loss  rate  of  the  flaming  zone.  In 
many  other  fires  there  was  considerable  glowing/ 
burning  char  left  after  passage  of  the  flaming  fire 
front.  Due  to  the  nature  of  our  fuels,  this  residual 
char  was  quickly  consumed  by  glowing  combustion, 
generally  leaving  no  residual  char.  This  trailing 
edge  of  glowing  char  obviously  had  very  little  effect 
on  the  propagation  of  the  fire  but  in  several  cases  it 
provided  a  significant  amount  of  the  mass  loss  and 
total  heat  released. 

So,  we  have  experimentally  measurable  reaction 
intensities  and  hence  by  equation  5,  the  reaction  ve- 
locity is: 

1    dw      .  -1  (6) 
p  _   —  .  — -  mm 
w  dt 

o 
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The  third  assumption  is  as  follows: 

3.  All  fuel  moisture  effects  on  reaction  velocity 
(and  intensity)  are  contained  in  a  separable  and 
independent  moisture  damping  coefficient,  r\M,  that 
is  defined  as  the  ratio  of  reaction  intensity  with  fuel 
moisture,  IR(Mf),  to  that  without  moisture,  IR(0). 
Each  experimental  fuel  bed  configuration  was 
burned  at  several  fuel  moisture  contents,  so  that  we 
have  experimental  measurements  of  r\M, 

IR(Mf) 

%  "    7^(0)  (7) 

to  which  we  try  to  fit  another  prediction  function, 

<r\u>=f(o,fr&Mf).  (8) 

Then  we  can  predict  the  value  of  the  reaction 
intensity, 

Moisture  Damping 

We  follow  Rothermel's  basic  concept  for  the  effect 
of  fuel  moisture  on  fire  propagation  and  for  the  defi- 
nition of  the  moisture  damping  coefficient  by  equa- 
tion 7  above:  that  is,  the  fuel  moisture  damping 
coefficient  is  the  proportional  reduction  of  reaction 
intensity,  relative  to  the  ovendry  fuel.  We  have, 
however,  revised  the  experimental  method  of  deter- 
mination and  reformulated  the  prediction  equations 
as  follows. 

Also,  the  experimental  determination  of  T|(M)  is 
actually  measured  as  a  ratio, 

R(Mf )  pbeQ(Mr )    R(Mf )  Q(Mf ) 

%  =  =  

%hvvroro  i?(0)Q(0)  (10) 

The  use  of  this  formulation  minimizes  measurement 
error,  explicitly  accounts  for  the  enthalpic  moisture 
load,  Q(Afp,  and  is  based  on  measuring  the  reduc- 
tion in  spread  rate  due  to  fuel  moisture.  (The  alter- 
native, used  by  Rothermel,  is  based  on  the  reduction 
of  mass  loss  rate  by  fuel  moisture  as  in  equations 
5  and  7.)  The  following  qualitative  observations 
are  noted: 

1.  In  most  of  the  excelsior  fires,  the  combustion 
zone  traveled  in  dry  fuel,  that  is,  except  for  the  very 
tightly  packed  fuel  beds,  the  fuel  was  dried  some 
distance  ahead  of  the  fire  front  and  no  free  water 
entered  the  combustion  zone.  In  these  cases,  the 
initial  fuel  moisture  content  affected  the  tempera- 
ture gradient  in  the  fuel  ahead  of  the  fire  and  hence 
affected  spread  rate  only  through  the  enthalpy, 
Q(Mp.  (This  is  similar  to  Albini's  conclusions  by 
computer  simulation  [private  communication].) 


Thus,  we  speculate  that  thermally  thin,  fine  fuels 
should  have  very  little  moisture  damping,  t\M  =  1, 
and  otherwise,  "moisture  damping"  may  imply  the 
ability  of  a  fuel  to  hold  and  carry  moisture  into  the 
combustion  zone  where  it  may  have  an  additional 
cooling  or  gaseous  dilution  effect. 

2.  A  second  reason  for  using  the  experimental 
method  of  equation  10  is  that  there  is  much  uncer- 
tainty and  ambiguity  in  defining  and  measuring  the 
"size  of  the  combustion  zone,"  which  area  is  divided 
into  the  mass  loss  rate  for  calculation  of  reaction 
intensity  (and  inherently,  the  mass  loss  rate  has  a 
larger  experimental  measurement  error  than  does 
spread  rate).  Also,  the  effect  of  the  enthalpic  mois- 
ture is  obscured  in  the  ratio  of  mass  loss  rate  over 
combustion  zone  area. 

If,  indeed,  some  fires  do  generally  travel  in  an 
"ovendry,"  moisture-free,  combustion  zone,  then  the 
reaction  intensity  must  be  independent  of  the  initial 
fuel  moisture,  that  is,  the  observed  effect  of  increas- 
ing fuel  moisture  would  be  to  slow  the  spread  rate 
and  decrease  the  size  of  the  combustion  zone  with- 
out, however,  changing  the  rate  of  mass  loss  per 
unit  area  within  that  zone. 

3.  In  the  above  mathematical  construction,  the 
two  methods  of  measuring  T|M  are  equivalent;  by  the 
spread  rate  method,  we  measure  a  more  dynamic 
variable  with  greater  precision  and  explicitly  ex- 
clude the  enthalpic  moisture  effect  from  moisture 
damping. 

In  this  new  formulation  of  moisture  damping,  the 
most  significant  improvement  is  to  completely  di- 
vorce moisture  damping  from  fire  extinction.  Each 
fuel  sample  must  still  have  a  characteristic  mois- 
ture content,  but  extinction  is  properly  handled  by 
the  probability  function,  P^.,  in  the  reaction  velocity, 
T,  discussed  later.  So  the  new  moisture  damping 
equation  is  very  much  simplified: 

<tim>  =  exp(-Mf/Mc)  (11) 

The  experimental  curves  are  shown  in  figure  la, 
a',  b,  c,  d,  e.  In  figure  la,  moisture  damping  (that 
moisture  effect  in  excess  of  Q(Mf)'s  effect)  is  shown 
to  be  negligibly  small  in  the  lightly  and  moderately 
packed  (P<0.04)  fine  fuels  in  support  of  the  observa- 
tions above.  Figure  la'  shows  the  moisture  damping 
(in  excess  of  Q{Mf)  for  the  very  compact  fuel  beds, 
|3>  0.04,  where  M.  =  0.50. 

Figure  lb  shows  a  definite  moisture-damping 
effect  in  Vi6-inch-thick  sticks  with  characteristic 
Mc  =  =  0.49.  This  relatively  large  Mj.  indicates  this 
fuel  is  easily  dried  and  thermally  thin.  Figures 
lc,  d,  and  e  show  the  progressively  greater  moisture 
damping  in  successively  larger  fuel  particles. 
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Figure  1a — Experimental  scatter  and  lack  of 
definitive  trend  for  moisture  damping  in  light- 
to-moderate  loads  of  fine  fuels  (excelsior). 


Figure  1a' — Moisture  damping  effect  in  the 
more  densely  packed — heavier  loaded  fine 
fuel  beds,  p>0.04.  The  curve  is  drawn  with 
characteristic  M.  =  0.50. 

c 


Figure  1b — Moisture  damping  for  a  =  25.2/cm 
(Vi6-inch  sticks).  Here  Mc  =  0.49,  indicating 
that  it  takes  a  considerable  amount  of  water  to 
dampen  fire  in  the  smaller  (thermally  thin)  fuels 
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Figures  1c,  d,  e — Moisture  damping  for  pro- 
gressively larger  fuel  particles:  a=  12.6,  6.3, 
and  3.15  cm-1,  respectively  (Vs-inch,  1/4-inch, 
and  1/2-inch).  The  characteristic  fuel  moistures 
are  as  follows:  Mc  =  0.23  for  a  =  12.6  cm-1, 
M=  0.19  for  a=  6.3  cm-1,  and  Mc=  0.15  for 
0  =  3.15  cm-1. 


Figure  2  shows  the  dependence  of  M,  on  fuel  par- 
ticle size  for  all  of  the  experimental  fuels. 

Equation  11  derives  from  a  dIR  =  -  (UMc)IRdMf 
type  of  model.  Note  that  the  characteristic  fuel 
moisture,  Mc,  is  shown  in  figure  2  as  functions  of 
only  the  size  of  fuel  particles;  the  experimental  fuels 
were  otherwise  very  similar,  all  being  sound  cured 
cellulosic  heartwoods — of  different  species,  however. 
An  empirical  equation,  Mc  =  /"(a),  could  be  con- 
structed from  the  data  in  figure  2,  but  would  be 
very  misleading.  Such  an  equation  would  only  sat- 
isfy pragmatic  needs  for  such  a  predictive  equation 
by  model  builders  of  fire  behavior  systems;  it  is  a 
strongly  held  opinion  of  the  author  that  values  of 
Mc  will  vary  greatly  with  the  physiology  of  different 
fuels  (bark  covered,  live  or  fresh-fallen  foliage, 


decayed  wood,  etc.).  This  opinion  is  based  primarily 
on  private  communications  with  Hal  Anderson  and 
extrapolating  his  moisture  response  observations 
(Anderson  1985)  into  a  forced  drying-fire  environ- 
ment. Hypothetically,  we  might  expect  that  a  fuel 
that  has  a  long  ambient  moisture  reponse  time  will 
also  more  easily  retain  moisture  by  forced  drying, 
thus  exhibiting  a  lower  Af  (smaller  values  of  Af,  will 
produce  a  stronger  response  in  T|M  due  to  change  in 
fuel  moisture  content,  Afp.  In  figure  3a,  b,  c,  we  use 
equations  10  and  11  and  replot  Anderson's  original 
data  from  which  Rothermel's  original  moisture 
damping  coefficient  was  derived.  Note  the  rela- 
tively strong  Af,'s  for  the  natural  pine  needles  when 
compared  to  cellulosic  pine  heartwood,  with  nearly 
the  same  surface  area-to-volume  ratio. 
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Figure  2— Plot  of  the  characteristic 
fuel  moistures,  A^,  vis-a-vis  fuel  particle 
size,  o.  Also  included  are  the  M^s  for 
Anderson's  pine  needle  data  from  figure  3. 
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Figures  3a,  b,  and  c— Moisture  damping  using 
Anderson's  pine  needle  data  and  "applying  the  new 
model.  The  characteristic  /Wc's  are:  Ponderosa  pine 
(3a)  Mc=  0.236,  lodgepole  pine  (3b)  Mc=  0.150, 
and  western  white  pine  (3c)  M  =  0.237. 
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Figure  4a — A  typical  curve  of  reaction 
velocity  as  a  function  of  the  packing  ratio, 
P.  This  plot  is  for  o  =  25.2  cm-1,  8  =  10 
cm,  and  pp=  350  kg/m3  pine  heartwood. 
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Figure  4b — Plot  of  calculated  reaction 
velocity  versus  measured  values  for  all 
experimental  fires.  Those  fires  in  the  up- 
per left  quadrant  were  experimental  fires 
that  did  burn  against  all  likelihood  and 
expectation. 


A  future  project  should  compare  Anderson's  am- 
bient moisture  response  times  with  these  M  's.  An 
attempt  should  be  made  to  relate  these  Af,'s  to  diffu- 
sion coefficients  and  other  physical  characteristics 
of  the  fuel  particles. 

Reaction  Velocity 

Note  that  Rothermel's  model  of  reaction  velocity 
was  characterized  by  an  optimum  packing  ratio 
which,  in  turn,  depended  on  the  fuel  particle  size; 
the  equations  predicted  maximum  reaction  velocity 
when  the  packing  ratio  was  at  its  "optimum"  value. 
Then  the  reaction  velocity  slowed  as  packing  ratio 
moved  away  from  optimum  in  either  direction,  but 
was  independent  of  fuel  bed  depth  and  hence  the 
total  amount  of  fuel.  The  reaction  velocity  as  de- 
fined here  has  a  similar  character  but  we  wish  to 
order  a  requiem  for  Pop.  In  the  extended  data  base, 
the  equivalent  effect  is  dependent  on  the  total 
amount  of  fuel  mass  and/or  fuel  surface  area  in  the 
fuel  bed.  Indeed,  as  the  fuel  bed  is  fluffed  up  to 
lighter  packing  (at  constant  load,  meaning  greater 
depth)  the  reaction  velocity  will  increase  as  fires  are 
shown  to  burn  faster  and  hotter  as  the  fuel  bed  is 
opened  up.  But  instead  of  reaching  a  maximum 
reaction  rate  and  then  slowing  down  as  a  full  con- 
tiguous fire  front,  the  fuel  becomes  sparse  and  the 
fire  front  breaks  up  and  ultimately  goes  out  as  de- 
termined by  the  function     (see  fig.  4a). 

The  fire  extinction  equation  that  provides  the 
practical  upper  limit  to  reaction  velocity  is  identical 
to  the  fire  state  (of  extinction)  probability  function 
in  Wilson  (1985): 

Pf  (nx)  =  1/(1  +  exp  (-7C  (nx  -  nx)  11.2  V~3)).  (12) 

Note  here,  however,  that  the  parameters,  1.2  and  nx 
(see  fig.  4a),  are  potential  "adjustment"  factors  that 
would  allow  T  to  work  correctly  for  a  different  fuel 
type;  a  number  smaller  than  1.2  will  "sharpen"  the 
edge  between  good  and  bad  burning,  and  conversely; 
nx  =  3.0  approximately  for  woody  cellulosic  fuels; 
larger  values  of  nx  will  move  the  "NO-BURN  cutoff' 
toward  heavier  fuel  loadings  or  lower  fuel  moisture. 
Note  that  nx,  given  in  equation  13,  is  strongly  de- 
pendent on  fuel  thermochemistry  and  fuel  moisture 
content. 

The  extinction  index  equation  that  partitions 
"good  burning"  from  "poor  burning"  fires  is  given  by: 

nx  =  ln(ai38  h0  IQw)KMf  +  Qf/QJ.  (13) 

Wilson's  extinction  index  is  described  in  detail  in 
Combustion  Science  and  Technology  (1985).  The 
index,  nx,  was  especially  designed  to  reflect  the 
energy  balance  of  fire  near  the  marginal  limits 
of  sustained  burning. 
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An  interesting  speculation  (without  experimental 
foundation)  is  to  include  the  wind  coefficient,  (J)^ 
in  the  argument  of  the  numerators'  logarithmic 
function;  in  addition  to  its  normal  effect  on  spread 
rate  (Rothermel  1972,  p.  26),  this  use  of  §w  also  has 
the  effect  of  "turning  r  on"  in  sufficient  wind  when 
otherwise  the  fuels  might  be  too  light  or  too  wet  to 
burn;  and  similarly  for  slope,  §s. 

We  must  caution  that  the  previous  formulation  of 
the  wind  coefficient  is  very  likely  to  produce  strange 
results  if  used  here  as  in  the  1972  fire  spread  model 
as  a  function  to  increase  rate  of  spread.  Cursory 
examination  of  the  "old  wind  model"  using  sketchy 
data  has  been  inconclusive  and  not  encouraging. 
Ralph  Nelson  (1988)  has  suggested  a  form  of  wind 
coefficient,  §  =  kw^'U^^lx  that  might  work  better. 

Until  a  new  wind  experiment  is  performed,  these 
formulations  should  be  used  with  great  caution. 
It  is  almost  imperative  that  a  new  set  of  laboratory 
wind  data  be  generated  for  revision  and  updating 
of  Rothermel's  wind  model  (and  similarly  for  slope). 
The  importance  in  the  field  of  wind  effects  and  the 
dearth  of  supporting  experimental  data  make  it  the 
next  priority  item  in  fire  behavior  research. 

The  new  experimental  formulation  for  reaction 
velocity,  r,  (fig.  4a  and  b),  is: 

<r>  =  0.34a  (apS)-172  exp  (-ap/3) 

Pf{nx)  min-1.  (14) 

The  0.34a  is  similar  to  Anderson's  old  8d  rule-of- 
thumb,  which  provides  that  the  burning  residence 
time,  t  in  minutes,  is  equal  to  8  times  the  fuel  par- 
ticle thickness,  d,  in  inches,  a  provides  the  Vd  de- 
pendence1 and  (a(38)  measures  the  absorption  cross 
section  of  fuel  particles  per  unit  area  of  the  combus- 
tion zone. 

The  exponential  function  of  ap73  limits  T  by  a 
characteristic  distance  that  propagating  flux  is  pro- 
jected through  the  fuel  bed  to  heat  fuel  ahead  of  the 
combustion  zone.  The  equation  components  with 
a(3  and  (38  arguments  slow  the  reaction  (velocity 
and  intensity)  for  heavier  and  denser  fuel  loads. 
The  Pf  function  shuts  the  fire  down  at  very  light 
loads  or  high  moisture,  for  example,  and  is  dis- 
cussed exhaustively  in  Wilson  (1985).  Fire  modelers 
in  the  future  will  probably  separate  this  probability 
function  from  the  reaction  velocity.  Only  in  very 
slow,  marginally  burning  fires,  with  very  deep  and 
very  lightly  packed  fuel  beds  where  flame  heights 
were  much  less  than  the  depth  of  the  fuel  bed,  did 


1The  new  data  can  be  made  to  fit  the  form,  (o"""1/p8)1/'' 
for  any  small  2<n<5;  it  seems  necessary  to  maintain  the 
1/d  dimensionality  and  provide  a  small  correction  to  the 
reaction  velocity  due  to  fuel  loading,  |}6  =  wo  /p  because 
only  the  fuels  near  the  surface  contribute  their  full  weight 
toward  propagating  the  fire,  particularly  in  deep  beds  of 
fine  fuel. 


the  Pf  function  contribute  significantly  to  the  burning 
rate.  Otherwise,  P^acts  only  to  turn  the  fire  on  or  off. 

Propagating  Flux  Ratio 

The  propagating  flux,  I  is  that  portion  of  the  fire's 
power  density  that  goes  to  drive  the  fire, 

Ip=Rpbe(Qf+MQJ.  (15) 

Following  Rothermel,  the  propagating  flux  must  also 
be  some  fraction,  h,,  of  the  reaction  intensity, 

IR,  such  that  ^  may  be  calculated  from  experi- 
mental measurements  by 

,       RphzQT    }  ■ 

hv  {dwjdt)  (16) 

and  an  empirical  prediction  equation  constructed 
to  fit: 

<^>=/,(a,p,5,--)  .  (17) 

It  is  only  fair  to  admit  here  that  all  of  the  parame- 
ters on  the  right  side  of  equation  16  are  measured  in 
the  experiment,  except  e.  Both  e  and  £  are  fractions 
that  diminish  their  respective  power  densities  so  that 
the  two  intensities  are  equivalent  to  one  another;  that 
is,  so  that  the  propagating  flux,  Ip,  equals  the  product 
of  the  supply  rate,  R,  and  the  heat  load  Q:  Ip  -  RpQ. 
It  follows  that  the  experimental  measurements  of  £ 
(effective  heating  ratio)  and  e  (effective  heat  number) 
are  interrelated.  As  a  consequence,  they  are  difficult 
to  measure  independently.  We  rely  on  Frandsen's 
(1973)  experiment  for  the  independent  prediction  of 
<£>  =  exp(-l/.22a).  Thus,  in  equation  17  we  are  actu- 
ally making  an  estimate  of  <£>/e  =  f  (a,p,8,-),  and  the 
ratio,  fye,  is  simply  a  rotation  operator  that  measures 
the  scalar  value  of  the  ratio  of  two  vectors  (the  hori- 
zontal propagating  flux  to  the  vertical  source  flux): 

8       hv  dwjdt    '  (18) 

In  any  case,  by  accepting  Frandsen's  expression  for 
e,  a  fairly  simple  straightforward  formula  for 
<£>  is  had. 

Note  that  we  have  here  an  interesting  experimental 
exercise  in  circuitous  reasoning: 

•  if,  in  equation  13,  we  interpret  R  =  dx/dt 

•  and  dwjdt  =  pb  '  d&dt 

•  and  assume  that 

dx  dj_  "flame  depth" 
d8      8  fuel  bed  depth 

then  the  experimental  measurement  of  IJz  becomes: 

.     _        df  _Pj£_ 

V  8  hv     '  (18a) 
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Figure  5 — Data  plot  of  £  versus  a(J. 


which,  factitiously,  compares  very  well  with  exper- 
imental determination  of  the  Zje  ratio  by  equation  18. 

Statistical  analysis  showed  that  the  product  of  the 
fuel  particles'  surface-to-volume  ratio  times  the  pack- 
ing ratio,  op,  is  by  far  the  most  dominant  variable 
governing  the  propagating  flux  ratio.  This  product 
is  dimensionally  equivalent  to  an  (optical)  absorption 
coefficient  of  power  transport  where  power  flux  den- 
sity decreases  with  distances  propagated  through  an 
absorbing  (heating)  medium.  This  formulation  is 
adaptable  to  a  model  of  propagating  energy  transport 
in  which  dlp  =  -klpdx,  and  k  =  cop,  etc.2 

Our  new  equation  for  propagating  flux  ratio,  <£> 
(fig.  5)  is: 

<£>  =  1  -  exp(-0.17  op).  (19) 

Note  that  the  factor,  0.17  cm,  (for  o  in  cm-1)  is  con- 
strued as  a  physical  skin  thickness  on  the  fuel  par- 
ticles that  is  absorbing  power  (being  pyrolyzed)  and 
is  comparable  to  Frandsen's  "effective  heating"  coef- 
ficient, 0.22  cm.  Susott  (private  communication)  has 
observed  a  charring-pyrolyzing  zone  in  woody  mate- 
rial that  consistently  measures  about  2  mm  deep. 

Furthermore,  we  should  note  that  all  of  our  experi- 
mental fuels  were  beds  of  long  narrow  particles, 


2The  new  data  will  also  fit  an  equation  with  fractional 
power  (that  is,  square  root)  of  cop ;  the  essential  character 
of  £  is  that  it  starts  at  zero,  then  increases  almost  linearly 
at  small  op,  then  rolls  off  at  large  oP;  the  theoreticians 
may  construct  a  theory  at  their  leisure. 


namely,  "one  dimensional,"  for  which  characteristic 
optical  absorption  distance  depends  on  the  geomet- 
ric cross-section  of  the  particle.  A  fuel  complex  of 
long  narrow  pine  needles  will  have  a  very  different 
optical  extinction  length  than  will  a  corresponding 
bed  of  wide  (two-dimensional)  oak  leaves.  The  point 
is  that  functions  that  depend  here  on  the  surface 
area-to-volume  ratio  of  the  fuel  particle  might  not 
be  universally  applicable  to  other  fuel  particle  con- 
figurations; particularly,  the  equations  of  propaga- 
ting flux  ratio,  reaction  velocity,  and  moisture 
damping. 

Rate  of  Spread 

The  predicted  value  of  spread  rate,  <R>,  is  given 
by  substitution  of  the  foregoing  equations  19,  14,  11, 
2  and  3  into  equation  1  (see  fig.  6): 

<>Kwo  <T><X\^ 

<R>=    •  (20) 

pbz(Qf  +  MfQw) 

From  figure  6  we  see  that  equation  20  overpre- 
dicts  rate  of  spread  in  the  faster  burning  fine  dry 
fuels  and,  to  a  lesser  degree,  in  the  very  slow,  light 
loads  of  larger  fuel  particles.  This  effect  is  the  re- 
sult of  combining  in  equation  20  the  several  pre- 
diction equations  for  <£>,  <r>,  <T1M>  which  were 
not  derived  from  statistically  independent  data 
sets.  Thus,  for  example,  the  same  experimental- 
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Figure  6 — Final  plot  of  the  distribution  of 
predicted  and  observed  rates  of  spread. 
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statistical  errors  in  spread  rate  are  applicable  to  <£>, 
<r>,  and  <i\i?>;  thus,  in  equation  20  they  are  applied 
thrice. 

Alternatively,  the  physical  problem  might  be  that 
the  particle-to-particle  variation  in  fuel  moisture 
or  heat  content,  etc.,  may  have  a  more  pronounced 
effect  on  spread  rate  in  the  drier  fuel  beds  with  few 
fuel  particles  (very  small  P).  These  variations  among 
the  individual  particles  were  not  measured. 

PROBABILITIES 

Let  any  one  of  the  critical  dependent  variables  such 
as  spread  rate,  R,  or  reaction  intensity,  IR,  be  repre- 
sented by  the  parameter  Y.  Then,  of  course,  the 
"real"  observed  value,  Y,  will  be  different  than  the 
predicted  fire  behavior,  <Y>,  which  implies  that 
there  is  a  probability  distribution  of  Y  about  <Y>. 

The  variable  <Y>  is  a  deterministic  function  of  the 
independent  variables,  (a,  p\  8,  Af,),  and  it  describes 
in  one  number  the  effects  (fire  behavior)  that  result 
from  the  conjunction  of  explicit  values  of  the  inde- 
pendent variables.  By  definition,  it  is  the  most  likely 
value  of  Y  that  may  occur,  but  it  is  very  unlikely  that 
Y  will  precisely  equal  <Y>.  It  follows  that  there  must 
exist  a  probability  distribution  of  the  real  fire  behav- 
ior when  we  are  given  the  expected  (predicted)  fire 
behavior  that  can  be  calculated  from  measured  val- 
ues of  a,  p,  5,  Mp  etc.  From  Wilson  (1987),  we  expect 
to  be  given  these  probability  distributions  so  that 
persons  that  must  take  risks  may  calculate  their 
chances  of  being  wrong. 

Now  let  us  review  the  probability  problem  and  so- 
lution. To  begin,  the  fire  manager  in  the  field  will 
have  at  hand  a  prediction,  <Y>,  of  predicted  fire  be- 
havior that  is  theoretically  different  from  the  actual 
"real"  behavior,  Y.  Also,  the  manager  will  have  some 
imperfect  knowledge,  P(Y),  of  the  prior  historical 
behavior  of  fire  in  similar  circumstances. 

The  manager  needs  to  know  the  probability  dis- 
tribution, P(Y/<Y>),  of  what  the  real  behavior,  Y, 
might  be,  given  prior  knowledge  of  Y,  P(Y),  and  the 
expected  behavior,  <Y>,  where  <Y>  is  estimated  from 
the  current  conditions,  a,  P,  6,       etc.  The  calcula- 
tion of  that  probability  follows  from  Bayes  theorem 
(Wilson  1987): 


P(Y/<Y>)  = 


P(<Y>/Y)  POT) 
P(<Y>) 


where  the  probability  distribution,  P(<Y>fY),  is  de- 
termined by  independent  experiments  such  as  the 
one  described  here. 

Figure  7  shows  the  error  distribution,  P(<R>/R) 
for  spread  rate,  <R>,  from  the  experiment  described 
here.  This  figure  is  the  one  to  be  used  for  probability 
calculations  of  rate  of  spread  while  using  the  predic- 
tion equations  presented  in  this  paper.  When  a  more 
precise  experiment  is  performed  or  more  accurate 
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Figure  7 — Plot  of  probability  that  a  rate  of 
spread  prediction  is  different  from  that  to  be 
observed.  Note  the  data  points  are  cumula- 
tive counts  and  fit  a  sigmoid  logistic  curve. 
The  bell-shaped  distribution  is  the  rescaled 
differential  of  the  cumulative  logistic  curve 
and  represents  the  expected  distribution  of 
the  predicted  rate  of  spread  about  the  true 
value. 


equations  are  derived,  then  a  new  P(<R>/R)  of  fig- 
ure 7  must  be  presented. 

These  error  distributions  give  a  clue  to  the  source 
of  the  "factor-of-two"  estimate  of  maximum  error  in 
fire  behavior  predictions  (Albini  1976);  we  note 
the  following  characteristics  of  these  distributions: 

1.  The  experimental  observations  of  <Y>/Y  are 
found  to  be  log-normal  for  all  of  the  dependent 
variables.  Thus,  observations  of  Y=  <Y>/n  and 
Y  =  n<Y>  occur  with  about  equal  probability. 

2.  The  standard  deviation  for  the  spread  rate 
observations  is  in  the  approximate  range  of 
(0.6<Y>)  2<Y<  (1-7<Y>).  (The  common  factor-of- 
two  error  would  be  about  two  standard  deviations 
on  the  log-normal  distribution.) 

3.  The  experimental  variation  of  the  independent 
variables,  a,  P,  8,  and  Mr  is  at  the  very  most  less 
than  5  percent  (rms  error  in  a  and  P  <1  percent; 

in  Mf  <2  percent  and  8  <3  percent. 

Thus,  the  experimental  errors  in  observation  of  the 
independent  variables  do  not  explain  all  of  the  vari- 
ance in  the  dependent  variables.  The  model  predic- 
tion error  is  most  likely  a  result  of  inaccuracy  in 
the  model  relationships  or  the  effect  of  missing  or 
neglecting  some  unknown  physiochemical -thermal 
parameter.  As  Albini  said,  we  "consider  models 
successful  if  the  relationships  predict  fire  behavior 
within  a  factor  of  two  or  three  over  a  range  of  two 
or  three  decades." 
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